Universal features of the order-parameter fluctuations : reversible and irreversible 

aggregation 
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We discuss the universal scaling laws of order parameter fluctuations in any system in which the 
second-order critical behaviour can be identified. These scaling laws can be derived rigorously 
for equilibrium systems when combined with the finite-size scaling analysis. The relation between 
order parameter, criticality and scaling law of fluctuations has been established and the connexion 
between the scaling function and the critical exponents has been found. We give examples in out- 
of-equilibrium aggregation models such as the Smoluchowski kinetic equations, or of at-equilibrium 
Ising and percolation models. 
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I. INTRODUCTION 

Fluctuations in many physical processes are difficult 
to analyze because they develop dynamically and often 
keep the memory of initial conditions. On the other hand, 
strong fluctuations are ubiquitous, as show examples of 
hadronization in strong interaction physics, polymeriza- 
tion, colloid aggregation, aerosol coalescence or the for- 
mation of large scale structures in the Universe. With the 
advent of recently developed advanced detection systems, 
the study of large fluctuations in physical observables be- 
came accessible in 'small systems', such as formed in ul- 
trarelativistic collisions of hadrons, leptons and nuclei, in 
the heavy-ion collisions at the intermediate energies or in 
the collisions of atomic aggregates [jl| . 

In theoretical studies, it is often assumed that fluctua- 
tions are irrelevant. In this spirit, many aggregation pro- 
cesses have been studied in the mean-field approximation 
. This problem has been revisited recently || . It was 
shown that contrary to the usual believe, fluctuations in 
the size distribution of largest cluster are generally large 
in the aggregation processes. Large fluctuations in the 
cluster multiplicity distribution have been also reported 
in the binary fragmentation process with the inactivation 
mechanism ||. 

This article deals with the features of fluctuations of 
physical quantities in a iV-body, <i-dimensional system, 
with N essentially finite. Moreover, the system is not 
necessarily at the thermodynamic equilibrium. Both 
these aspects of our approach are important in many ar- 
eas of physics where, for example, small and strongly 
fluctuating systems are produced in the violent collision 
processes. Consequently, these systems live shortly and 
the typical time-scales are such that standard methods 
of equilibrium statistical physics may not to applicable. 

We shall be particularly interested in self-similar sys- 



tems, such as the fractal objects or thermodynamic sys- 
tems at the second-order phase transition. The self- 
similarity means in particular that one is unable to de- 
fine the characteristic length ~ (N*) 1 ^, where N* is the 
characteristic size, which could be associated with the 
disappearance of fluctuations. Our aim in this work is 
the discussion of universal scaling laws of fluctuations of 
different observables in self-similar systems. In particu- 
lar, we shall consider the order parameter fluctuations in 
any system, both equilibrium and non-equilibrium one, in 
which the second-order critical behaviour can be identi- 
fied. These considerations will provide an understanding 
of the relation between the order parameter, the criti- 
cality and the scaling law of fluctuations. The notion of 
relevant variable (relevant observable) for the discussion 
of critical behaviour in finite systems will appear from 
this discussion. 

The paper is organized as follows. In Sect. II the or- 
der parameter fluctuations in statistical systems are ana- 
lyzed and the relation with the finite-size scaling analysis 
and Widom's hypothesis is developed in details [[| (Sect. 
II. B). The generic features of the tail of the scaling func- 
tion is addressed in Sect. II. C. In Sects. II. E - II. H, 
the generalized scaling of the observable quantities , the 
A-scaling, is discussed and the reasons for the deviations 
from the limiting cases A = 1 and A = 1/2 are presented. 

Sects. Ill - VI are devoted to the detailed discussion of 
several well known generic models, using the results and 
methods of analysis proposed in Sect. II. In Sect. Ill, 
the non-critical fragmentation model is discussed which 
exhibits the power-law cluster size distribution. Results 
obtained in Potts model are discussed in Sect. IV. Sect. 
V is devoted to the discussion of the reversible aggrega- 
tion, as modelled by the percolation model. Both real- 
istic 3-dimensional bond percolation and the mean-field 
percolation on the Bethe lattice is being considered. Ir- 
reversible aggregation and the fluctuation properties in 
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the Smoluchowski kinetic model are discussed in Sect. 
VI. Finally, the main conclusions are given in Sect. VII. 

II. ORDER PARAMETER FLUCTUATIONS 

A. Correlation function argument 

Let us call m the observable under investigation. For 
the reason of presentation, we shall restrict ourselves 
to the case where to is a scalar quantity and takes 
real positive values. ( In fact, this restriction is not a 
true limitation and, generally, one can consider \m\ as 
well.) Fluctuations of the order parameter in thermo- 
dynamic systems are expected to have different prop- 
erties at the critical point and outside of it. Far from 
the critical behaviour, the correlations are short-ranged. 
Fluctuations of the extensive order parameter to in this 
case resemble the ergodic Brownian motion of this vari- 
able in its proper configuration space. Consequently 
< (to— < to >) 2 >/< m > is roughly a constant, mean- 
ing independent of the number of constituents in the 
sample. On the contrary, close to the second-order tran- 
sition point, the fluctuations are correlated throughout 
the whole system and the correlation length £ for the in- 
finite system becomes infinite as well. Let us now define 
in such a system the deviation e of the driven parameter 
from its critical value, to the local value of the order pa- 
rameter and the field h conjugated to to. The following 
isotropic correlation function is introduced : 

a n (e,h, ri,..., r n -i) = < m(r a )m(r a - f\) . . . 

■ ■■m(r - r n -x) > . (1) 

Notation < ... > in the above expression denotes the 
thermodynamic average at a given e and over the position 
r . For vanishing e, all the length scales disappear and 
the correlation length £ diverges algebraically as 

with a universal exponent v, which depends only on the 
universality class of the transition. Moreover, the scal- 
ing description of second-order critical phenomena leads 
to the fundamental postulate that the thermodynamic 
potential G verifies 



G(Ae, A 



2-CH-/3 



h) ~ X 2 - a G{e,h) 



(2) 



and this defines the two universal exponents a and [3. Let 
us now go back to the correlation function a n defined in 
(|l|). The integral of a n over the n — 1 space- variables 
: fi, ...,r„_i, is equal to the n-th derivative of G with 
respect to the field h. Hence, to be consistent with (||), 
we must get the scaling relation : 

a n (Ae, \ 2 - a -?h, \~ v f\ , . . . , A-V„_!) ~ 

v n/3 



For n = 1, we recover the well-known scaling behaviour of 
the averaged order parameter with the critical exponent 
p. For any integer n, putting all the space-variables to 
in formula (||) leads to the scaling of powers of the local 
order parameter : 



< m™ > (Ae, X 2 - a - p h) ~ A n/3 < to" > (e, h) 



(4) 



As a consequence, if we let h = and A = 1/e, one finds 
that the quantities < to™ > / < to >™ are independent 
of e, when close to the transition. This means, by the 
finite-size analysis, that this ratio is independent of the 
size N of the system at the transition point. 

Let us now introduce the cumulants n q from the gen- 
eral formula for the moment expansion of the order- 
parameter probability distribution P[m] : 



In [ P[m] exp(mu) J = — -n q 
\m=0 / g=o 



(5) 



Expanding the l.h.s. of the above expression in the power 
series in u, and comparing corresponding powers on l.h.s. 
and r.h.s., one derives the relations between ordinary mo- 
ments of P[m] and the cumulant moments : 



Ki = < to > 

K 2 = < TO 2 > — < TO > 2 

K3 = < TO 3 > -3 < TO 2 >< TO > +2 < TO > 3 

K4 = < m A > —4 < to 3 >< TO > —3 < to 2 > 2 

+ 12 < TO 2 >< TO > 2 -6 < TO > 4 



(6) 



In the case of the second-order phase transition, as a 
result of scaling relations (^) , all cumulant moments scale 
like 



Consequently, the generating function of the m- 
probability distribution (|^) is only a function of the re- 
duced variable < to > u. This can be written as : 

1 f 27r - 

-F"[ TO ] = 7T~ / G(iu) exp(—imu)du (7) 
2tt Jq 

where G is the generating function : 

oo 

G(u) — P[m] exp(TOu) . 

m=0 

In the case when < to > tends to oo but to/ < to > 
remains finite, we can rewrite formula (Q) in the more 
compact form : 



A™^<r n (e, /i,r*i, ...,r„_i) . (3) 



< to > P\m] = 



[m] = * (-^-) 

V< TO >/ 



(8) 
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which is valid at the critical point. <E> is the scaling func- 
tion of the single reduced variable mj < m >. As stated 
before, we can express this scaling as : 



E 



P/v[to] exp(mw) = W(< m > u) 



(9) 



which is the necessary and sufficient condition for the 
applicability of scaling law m) [M . This result implies 
also that if this scaling occurs for the fluctuations of the 
parameter m then it holds also for fluctuations of any 
power : 

X = N a m b 

of this parameter as well. This is a consequence of (^) 
and 

P N [X]dX = P N [m]dm . 

Till now, we did not specify reasons of changing < 
to >. Indeed, under the condition that the scaling frame- 
work of second-order phase transition holds, the scaling 
relation (||) is valid independently of the explicit rea- 
sons of changing < to >, and independently of any phe- 
nomenological details. In other words, the explicit rela- 
tion between the size N of the system and < to > need 
not to be known at this stage. In the following section, 
we shall show how to derive supplementary informations 
about $ when the system is at the pseudo-critical point. 



B. Widom's hypothesis and the finite-size scaling 
argument 

The hypothesis of Widom states that in the thermody- 
namic limit of a system at thermal equilibrium, the free 
energy density close to the critical point scales as M : 



/(A^,Ae)~ X 2 - a f( V ,e) 



(10) 



where a,/3 are the usual critical exponents, 77 is the in- 
tensive order parameter and A is the scale parameter. 
Even though finite systems do not exhibit the critical be- 
haviour, nevertheless their properties may resemble those 
of infinite systems if the correlation length £ is larger or 
comparable to the typical length L of the system. In 
this case, one speaks about the pseudocritical point in a 
finite system at a distance 



cN- 1/vd 



(11) 



from the true critical point [||. The quantity N in ( Jlj] ) 
is the size of the (i-dimensional system and c is some 
dimensionless constant which can be either positive or 
negative. This constant c is negative if the maximum of 
finite-size susceptibility or of any other divergent macro- 
scopic quantity in the thermodynamic limit lies in the 



ordered phase, while c is positive if this maximum is in 
the disordered phase. One can then derive the finite-size 
scaling of the total free energy : 

F(r,,e,N)=Nf(r 1 ,e) 

at the pseudocritical point : 

F c (r),N)~f(r,N&) . (12) 
In deriving ( |l2"| ) we used the hyperscaling relation : 
2 — a = vd . 

The canonical probability density of the order parameter 
p n[v] is g iven by |§ : 

PN[ri\ = ^-exp(-0TF(ji,e,N)) , (13) 

where the coefficient /3t(= l/T) is independent of 7/ (T is 
the temperature of the system) . Using Eq. ( p"3| ) , one may 
calculate the most probable value of the order parameter, 
which is the solution to the equation 



dP N [y] 

drj 



= 



as well as the average value of the order parameter and 
the partition function 



Z 



N 



N~ — 



<rj> 



(14) 



rf in (14) denotes the most probable value of the order 
parameter. The average value of the order parameter 
vanishes for large values of N, since both (3 and 2 — a — 
2/5 + 7 are positive. The probability density Pn M obeys 
then the scaling law formally identical to Eq. IM) : 



< 77 > PnItj] 
where, in addition, 



V 



<T]> 



$(z) - exp (-(3 T f(az,c)) 



(15) 



(16) 



In the above formula, we have omitted the temperature- 
dependent multiplicative factor which can be determined 
by normalization of Pn [t]] ■ Coefficients a and c may both 
depend on j3r ■ We can then rewrite the scaling ( |l5| ) in a 
standard form for the extensive order parameter m = Nr/ 



< m > P N [m] = $(z(i)) 
with the scaling variable Zm defined by 



Z(l) 



TO — TO 



< TO > 



(17) 



(18) 
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m* denotes the most probable value of the extensive or- 
der parameter. We call ( [lT| ) with (|l8|) the first-scaling 
law. The scaling domain is defined by this asymptotic 
behaviour of Pjy[m] when to — > oo and < m >— > oo, but 
Z(ys has a finite value. The normalization of the proba- 
bility distribution P/v [v] an d the definition of the average 
value of to, provide the two constraints : 

/oo 
§{z (1) )dz {1) = 1 
-m* /<m> 

lira / z (1 )<i>(z( 1 ))dz( 1) = 

-m* /<m> 



The first-scaling law (15) is a consequence of the self- 
similarity of the statistical system. The self-similarity 
means here that the fluctuations of the reduced order pa- 
rameter r\j < i] > at different scales characterized by 
different values of the intensive order parameter < rj >, 
have identical properties. This is a qualitative explana- 
tion for this scaling. 

The logarithm of scaling function ( pf ) corresponds to 
the non-critical free energy density at the renormalized 
distance e = c from the critical point. If it happens that 
the order parameter is related to the number of frag- 
ments, like in the Fragmentation - Inactivation - Binary 
(FIB) process piO , then ( |l5| ) can be written in an 
equivalent form to the Koba - Nielsen - Olesen (KNO) 
scaling , proposed some time ago as an ultimate sym- 
metry of the S - matrix in the relativistic field theory 
||l9f . The multiplicity distribution of produced parti- 
cles is intensely studied in the strong interaction physics 
where simple behaviour of much of the data on hadron- 
multiplicity distribution seems to point to some univer- 
sality independent of the particular dynamical process. 

If, instead of real positive scalar, the parameter un- 
der investigation is a vector of dimensionality n : to = 
[mi, • • • , TO n ] , then the first-scaling law ( |Tt| , |l8| ) takes the 
more general form : 



<\m\ > n P N [m] = 



with 



z (i) 



TO 



TO 



< TO > 



(19) 



(20) 



The scaling limit in ( |l9| ) is defined by the asymptotic 
behaviour of P/v[m] when m, — > 00 (i = 1, • • • , N) and 
< \rh\ > — > 00, but zu\\ (i = l,---,N) have finite 
values. 



C. The tail of the scaling function 

The scaling function $ introduced in the first-scaling 
law ( fT^ ) has some typical features reminiscent of the non- 
Gaussian critical distribution of the order parameter. In 



this section, we are interested in the behaviour of scaling 
function for large values of the reduced parameter to/ < 
to >, so we have to study the system subject to the 
small field h conjugated to the order parameter. This 
breaks the symmetry of the distribution by shifting to 
and < to > towards larger values. More precisely, let us 
write the probability to get the value 77 of intensive order 
parameter at the distance e from the critical point as : 



P N [r], e, h] = Pn[v, e > °] exp((3 T r]Nh) 



(21) 



Till now, we have studied the behaviour of Pjy[r),e,0] 
for which the first-scaling law holds when e = (the 
critical point) or £ = cN~ x l vd (the pseudo-critical point). 
Substituting ( |l5| ) and using (|lj), we obtain : 

P N [r), 0, h] = exp (ln$ (yN^^j + /3 T rjNh^ . 

The most probable value 

of the order parameter in the limit of small external field 
h, is given by the maximum of the term in the exponen- 
tial. Since 

. 2-a-/3 
S= — 3 — 

we get : 

ln$(/i7iV^r) (3 T 0Nh . (22) 

Relation ( p2| ) is valid for any value of N if and only if 

$(2) - exp(-az 5+1 ) = exp(-a/) , (23) 

with the coefficient a which depends on the temperature 
regularly. 

One can express this relation in a different way. The 
anomalous dimension for extensive quantity to = Nr] can 
be defined as : 



q = lim qm = lim 

at^oo^ jv-»oodln7V 



(In < to >) 



(24) 



One can see that due to both ( [14[ ) and the Rushbrooke 
relation between critical exponents : 

a + 2/3 + 7 = 2 , 

the anomalous dimension is : 





1 



7 + 2/3 



(25) 



Since both a and [3 are positive, therefore g is contained 
between 1/2 and 1 for equilibrium systems at the criti- 
cal point of the second-order phase transition . Because 
of these additional relations between critical exponents, 
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one may note that a behaviour of the tail of the scaling 
function (B3h is governed by the exponent : 



1 



1 



= 6 + 1 = 



(26) 



which is always larger than 2. The limiting case : v = 2 
, i.e. the Gaussian tail (see (|23|)), corresponds in this 
framework to the non-critical system. 

Finally, let us mention in passing that whenever 
cluster-size can be defined in a system exhibiting the 
second-order phase transition , like e.g. in the case of 
percolation, Ising model or Fisher droplet model, the ex- 
ponent t of the power-law cluster-size distribution 



satisfies additional relations 



7 + /3=- 
a 



7 + 2/3 = 



what yields : 



and 



r - 1 



(27) 



(28) 



Since g at the second-order equilibrium phase transition 
is contained between 1/2 and 1, therefore the allowed 
values of exponent r at the critical point are : 2 < r < 3, 
and the normalized cluster-size distribution : 

JY 



^2 sn s = N 



s=l 



Ns- 



Consequently, whenever r is defined, we get an infor- 
mation whether the studied equilibrium system is at the 
second-order phase transition and whether the consid- 
ered extensive quantity can be identified with the order 
parameter of this transition. 

Let us define, for example, the multiplicity as the total 
number of clusters. (The definition of clusters includes 
also the monomers.) The cluster multiplicity cannot be 
an order parameter of these kind of equilibrium phase 
transitions because with 2 < r < 3 : 



N 



N 



what means that the average multiplicity scales as the to- 
tal mass of the system at the transition point, i.e. g = 1. 
On the other hand, the size of the largest cluster is a nat- 
ural order parameter for these kind of phase transitions. 
In this case we have : 



< Sr, 



(29) 



what is a direct consequence of : 



N 



n * ~ 1 



s=<s max > 

i.e. that there is in the average only one largest cluster. 
Moreover, the relation ( p7[ ) derived for the second-order 
critical phenomenon is correctly recovered. One should 
emphasize, that the relation (^9|) is very general and its 
derivation does not depend on the assumption of thermo- 
dynamic equilibrium. In other words, the relation (|29| ) 
between the anomalous dimension and the exponent r 
, is valid also for the off-equilibrium second-order phase 
transitions. We shall return to this point in Sect. VI. 

The cluster multiplicity could be the order parameter 
whenever r < 2, though this cannot happen in the equi- 
librium phase transitions. Note that this argument to 
sort among different candidates for the order parameter 
, requires only the knowledge of r, i.e. the complete in- 
formation about the critical process is superfluous. We 
shall use this argument later in the case of percolation 
model and Smoluchowski model of gelation. Finally, we 
shall see below in the Mekjian model that we may have 
power-law size distribution with r < 2, in the absence 
of phase transition governed by the multiplicity as the 
order parameter. 



D. Landau - Ginzburg theory of phase transitions 

Let us consider the Landau - Ginzburg (LG) theory 
as an exactly solvable example of the second-order phase 
transition. The homogeneous LG free energy density is : 

where b is a positive constant. The most probable value 
of the order parameter rj in the disordered phase (e < 0) 
is implicitly set to 0. It is more convenient to work with 
the extensive order parameter m = Nrj when dealing 
with the finite systems. The probability of a state m for 
a given e is |j) : 



P N [m] = — exp 
am 



? 4 

, m ,m s 



(30) 



Zn is defined by the normalization of Pm [m] ■ To remain 
consistent with other sections of this paper and without 
loss of generality, we consider now the case where m is 
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positive. We will admit that A is so large that the first 
two terms in the free energy expansion are sufficient to 
study the phase transition. At the critical point : e = 0, 
the leading term of the free energy density is propor- 
tional to to 4 . Standard integrations yield the values for 
the partition function Zjq and the average value of the 
order parameter < to >, both proportional to TV 3 / 4 . In- 
troducing them in (J30|), one finds : 



x exp 



T 4 [l/4] 



in 



< to > 



(31) 



which has the form of (g) . Note that the scaling function : 
$(z) exp(— z 4 ), decreases very fast as one moves away 
from the most probable value. This result is consistent 
with the analysis done in the previous section. 

The pseudo-critical point is the value of e for which the 
finite-size thermal susceptibility reaches its maximum. 
Writing that the inverse of this susceptibility is the sec- 
ond derivative of the free energy with respect to the order 
parameter, one finds : 



e = -6 



£[3/4] 
r[l/4] \f3 T N 



1/2 



(32) 



This results is correct at the first order in AT x l 2 . Replac- 
ing e in (p0[) by (p2j), leads to the scaling form of P/v[to] 



< to > Pjy [m] = A exp 



r[i/4] 

6(m/ < to >) 2 )] 



r[3/4] 2 , , , 4 

L ' 2 ((TO/ < TO >) 4 



(33) 



where A denotes a normalization constant. We recover 
indeed the first-scaling law with the exponential tail : 
exp(— az 4 ), for the large arguments. 

Outside of the critical point in the disordered phase 
(e > 0), the leading term of the free energy is propor- 
tional to m 2 , and the probability distribution Pjv[to] is 
essentially Gaussian. Deriving, as previously, the values 
of Zjy and < to > (both behave like A 1 / 2 in this case), 
we get the scaling form : 



r l 4 f 4 

< ?7i > Pn [to] — — exp 

7T V 7T V < TO > 



(34) 



which is still under the form ( |l5| ) but with a Gaussian 
scaling function reminiscent of the Gaussian fluctuations. 

Finally, in the low temperature regime (e < 0), the 
most probable value of the order parameter is positive : 



-3* 



Developing P/v [to] in ( |30| ) around this point leads to the 
expression : 



in " 4 / 2 P/v[to 



x exp e* — 2 



} e (to — m*) 2 
b to* 



(35) 



which is no more in the standard form (^). In this case, 
the average value of the order parameter < to > is of the 
same order of magnitude as its most probable value to* 
and one can rewrite (pq) in the scaling form : 



< to > 1/2 P/v[to] ~ exp ( -a 



(to — TO*)' 
< TO > 



(36) 



where a is a positive constant. This particular scaling 
form will be discussed later in details. 



E. The A-scaling law 

One may ask, what happens if the observable quantity 
is not the order parameter but the A-dependent function 
of the order parameter like : 



where 



TO = A" 1 - A a2 77 



ax > g + a% — 1 



(37) 



(38) 



The latter condition ensures that the order parameter 
does not determine the leading behaviour of to. For large 
A : 

< to > - A Q1 . 

Writing ( |i"7| ) with m instead of 77 and taking into account 
that : 

Pjs[[rj\dri — P/v[to]g?to 
one finds the generalized law : 

< to > A P N [m] = <fr(z (A) ) = $ ( m ~ m * A ) , (39) 



< TO > 



where : 



A= g + a2 ' 1 <l . 
ax 

This generalized law will be called in the following the A- 
scaling law. The scaling function <f>(zrM) depends only 
on one scaled variable : 



S(A) 



< TO > 



(40) 



The normalization of the probability distribution P/v[to] 
and the definition of the average value of m provide two 
constraints : 



G 



lim 



lim 



$(z A )dz A = 1 



za§ {z A )dz& = 



which are consistent with : A < 1, because the scaling 
function $ is positive. The scaling function $(z(A)) i n 
( |39| ) has an identical form as except for the inver- 

sion of the abscissa axis. In particular, its tail for large 



2(A) 



has the same form : 



$0 (A) ) ~ exp (-*f A) ) - cxp (-zffij (41) 

as given in (p3[). One should mention in passing that if 
ai < g + <Z2 — 1 



in (37), then 



< m >~ Af a2 < 77 > 

and A = 1, following the remark of Sect. II. A. 

According to (|l6|), the logarithm of scaling function 
$0(A)) : 

ln$(z (A) ) = -#r/(a*( A ),c) 

is related to the non- critical free energy /, in either or- 
dered (c > 0) or disordered (c < 0) phase. 

As an important example, we see from (37), (|39|), that 
the A-scaling of the extensive variable : 

to = N(l - 77) = Nfj 

can be used to determine the anomalous dimension since 
in this case : A = g. For this reason , to is a very useful 
variable in all phenomenological studies. At the phase 
transition : 

<Nr}> ~ N 
but the finite-size corrections are algebraic. 



F. Off-critical scaling 

A — 1/2 with $(z(a)) nearly Gaussian, is a particu- 
lar case of a A-scaling associated with the non-critical 
systems This limit : 



< m > 1/2 -Pat [to] = <f> 



< to > : / 2 
*(*(i/a)) - 



(42) 



which is called the second-scaling law, has been found in 
the shattering phase of the non-equilibrium FIB process 



pq| and in the 'gaseous' phase of equilibrium percola- 
tion process ^3|. We should recall also that this form of 
scaling function has been seen for LG model in the low 
temperature regime (see (|36|)). 

More generally, let us suppose now that the extensive 
parameter to is not critical, i.e. either the system is in 
a critical state but the parameter to is not critical, or 
the system is outside of the critical region. The value of 
to at the equilibrium is obtained by minimizing the free 
energy. The free energy F is analytical in the variable to 
close to its most probable value to* : 

F~ Ar _1 (m-m*) 2 . (43) 

Using ( f43f ) one obtains 

< to >~ li*N , 

where fx* is a positive (finite) number independent of A^ 
and : 



N l/2 



< TO > 



1/2 



(44) 



The probability density Pjsr [to] verifies the second-scaling 
law (p9|) : 



< to > 1/2 [to] = exp |^-/3 T M* 

= *(*(l/2)) ■ 



TO— < TO > 
< TO >!/2 



(45) 



This is a particular case of the A-scaling law (A = 1/2) 
and the scaling function is now Gaussian [2^] . This scal- 
ing ( f42| ) holds for < to >~ A^ but now with the exponen- 
tial finite-size corrections. This is a principle difference 
from the finite-size corrections and/or A-scaling. The 
above arguments apply to any second order phase tran- 
sition. In particular, they are not limited to LG theory 
of phase transitions (see Eq. (|36|) 



G. Finite-size cross-over effects 

The discussion of previous section is valid for systems 
at the critical (and pseudo-critical) point or far from the 
critical point in the ordered phase. Let us suppose now 
that the system is prepared such that : 



< to >~ N 9 



9 <1 



and g is not the anomalous exponent. Here, we would 
like to study how the finite system evolves when the con- 
trol parameter e tends slowly to 0, namely : 



N 2g 



We shall address this question in the mean-field approx- 
imation using LG theory. Let us first write down the 
average value : 
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< to >= 



J °° mexp(—em 2 /N — bm 4 /N 3 )dm 
J °° exp(— em 2 /N — bm 4 / 'N 3 )dm 



(46) 



Hence, writing this definition with the new driving pa- 
rameter e' = eN 1 / 2 and using the rescaled variable 
ml = m/N 3 / 4 , the average value of m can be put in 
the form : 

< m >= N^^eN 1 ' 2 ) , 
while its most probable value is : 



If the exponent g is not too small, i.e. if e does not 
vanish too fast, the two quantities : < m > and m*, 
have to coincide. This is because the exponential weight 
term in (Efj) diverges as ~ exp(e 2 iV/4&) when 



eN 1/2 ~ N- 



(47) 



becomes large with increasing TV. As a consequence, the 
common behaviour of < to > and m* is : ~ N 9 . The 
scaling form (fiff) in this case is : 



< to > A Pjsi[m] ~ exp I — c 



(to — m*y 



< m > 



2A 



(48) 



with c a positive constant, and : 

2 y 

We recover here the two cases previously discussed in 
Sect. II. B. When g = 1, i.e. when e = const, then this 
is the second-scaling law. When g =3/4, then this is 
the first-scaling law since the finite system is yet in the 
critical region (g = g). In between these two limiting 
cases, A-scaling holds with 1/2 < A < 1. Note also that 
the scaling function in (^) has a Gaussian form, even 
for A > 1/2, which is quite different from the case ( |2^ ) 
of Sect. II.C. 



H. Summary : panorama of the A-scaling for 
thermodynamic systems 

Several features of finite systems are important if one 
wants to study either the criticality of the correspond- 
ing infinite system or the distance to the critical point. 
One should name here : the A-scaling (this includes the 
first-scaling law A = 1 as well), the form of the tail 
of the scaling function <E>, and the anomalous exponent. 
All these features are closely related with the properties 
of the scaling function which characterizes finite system 
at the equilibrium. If the infinite system experiences a 



second-order phase transition, and if m is the scalar or- 
der parameter or the shifted scalar order parameter (|3^) , 
then : 

At the critical point, the corresponding finite system 
exhibits first-scaling law if m is the order parameter, or 
A-scaling law if m is the shifted order parameter. In 
both cases, the tail of the scaling function ~ exp(— z v ) 
is characterized by a large value of the exponent v — 
1/(1 — g) > 2, with g being the anomalous exponent, 
i.e. the exponent characterizing decrease of the extensive 
order parameter with the size N of the finite system. 
The values of A are restricted to < A < 1, and the 
anomalous exponent g takes values in between 1/2 and 1 
for second-order at-equilibrium phase transition. 

Far from the critical point, finite system exhibits 
second-scaling law with the Gaussian tail of the scaling 
function. 

Close to the critical point when e — > if N — + oo, the 
finite system exhibits the cross-over phenomenon from 
the first-scaling to the second-scaling law by a contin- 
uous A-scaling law with Gaussian shape of the scaling 
function. One should remind here that the precise depen- 
dence of e(N) is irrelevant provided that : g < g < 1, 
i.e. that the conditional point approaches not faster 
than the pseudocritical point. This last remark is impor- 
tant in the phenomenological applications of the scaling 
theory to the situations where the N - dependence of 
the conditional point is governed by an external control 
parameter with unknown relation to the system size. 

Last but not least, if the parameter to is not singular 
at the transition, then all properties of its probability 
distribution are the same as in the case of non-critical 
systems. 

In phenomenological applications, it is often difficult to 
get the probability distribution with sufficient accuracy 
for the values of scaling variable which are far from the 
most probable value since this corresponds to very small 
probabilities. It is then more judicious to work with mo- 
ments of the distribution instead of the distribution itself. 
For example, when the system undergoes A-scaling, the 
properly normalized cumulant moments (||) : 



(ki)« a 



const 



(49) 



are independent of the size of the system. An important 
consequence is that the generating function of the re- 
distribution : G(u) — ^2 P[m] exjp(rnu) , is a function 
of the reduced variable < m > A u only, generalizing a 
remark of Sect. II. A for the generating function in the 
first-scaling case. 
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III. A NON-CRITICAL MODEL : THE MEKJIAN 
MODEL 

The Mekjian fragmentation model is an equilibrium 
model which describes the decomposition of system into 
an ensemble of fragments. The statistical weights for 
every configuration of fragments are given explicitly in 
this model. If n s denotes the number of fragments of size 
s with the size-conservation : JV — J2 S sn S' the weight 
function for the configuration {n s } is given by fll3[] : 



N 



W N ({n s },x) = ]J 



L n s ls ns (x + s — 1) 

with x being a real control parameter. Many exact re- 
sults can be obtained in this simple model. Here, we are 
interested in the multiplicity distribution Pjv[m], where 
the fragment multiplicity is : m = n s . We can show 
that 111 : 



P N [m]=x m \S 



(m)| 
JV 



T(x) 



T(N + x) 



where |S , jv™' ) | are the signless Stirling numbers of the first 
kind. Knowing then the generating function for these 
Stirling numbers : 



l in e ' 



m=0 



T{x)T(xe u + JV) 
T(xe u )T(x + JV) 

one obtains the average value of to : 

1 



(50) 



< to > 



N 



xlriN + (x- 1)7- tp(x) +0(1/N) 



Moreover, making an asymptotic development of (50) for 
large JV and small s, one obtains : 



rn— 



P N [m]e mu ~ N x{eU -^ . 



(51) 



The latter approximation is known to be correct for finite 
values of u Jl5"[ . This means that P/v [to] is approximately 
a Poisson m-distribution with parameter x In JV. In the 
leading order we have then : < m > ~ to* . Inverting 
( pp| ) to get Pjv [m] as a Fourier transform, and making JV 
large yields the scaling formula : 



< to > 



exp 



- (x - 1)(7 - ip(x)) 

o- 1 



(m— < TO >) 2 
2 < TO > 

TO— < TO > 



< TO > 



< TO > 



(52) 



This is nothing else but the second-scaling law Q42 
for the multiplicity distribution when JV becomes large 
enough, because < to >~ m* . When < to > is large 



enough, the second term in (52) is always very small com- 



pared to the first one for a finite x. 

Different fixed values of the control parameter x mimic 
different situations of the fragmentation. For x -C 1, one 
has the situation of a fused system. For x ~ 0.5, the frag- 
mentation resembles the evaporation of light fragments. 
The limit x 1 corresponds to the complete dissocia- 
tion of the mass into light fragments (monomers). Each 
of this situation is characterized by a different fragment- 
size distribution. The case x = 1 is particular in this 
model since it leads to the power-law size-distribution 
( p7| ) with the exponent r = 1. Following discussion in 
Sect. II. C, cluster multiplicity could be the order pa- 
rameter. On the other hand, the second-order equilib- 
rium phase transition is associated with 2 < r < 3, what 
implies that the equilibrium model of Mekjian is a non- 
critical model. Indeed, that is what can be seen also 
in the cluster- multiplicity scaling law (|52"|). Hence, the 
power-law cluster-size distribution alone does not guar- 
antee that the system exhibits the critical behaviour of 
any kind |0J. 



IV. EXAMPLE : THE POTTS MODEL 

A generalization of the magnetic Ising spin model has 
been proposed by Domb [Q and studied in details by 
Potts jl(|. In this model, one considers a system of JV 
sites in the d-dimensional space. The magnetic state of 
each site i is characterized by a quantity called : a spin 
(say : s,-). Each spin is of the same constant modulus and 
points to one of the q equally spaced directions, labelled 
from to q — 1. The ferromagnetic short-ranged Potts 
Hamiltonian is then : 



H, 



(53) 



where 5 is the Kronecker symbol, and J is the posi- 
tive coupling constant. The sum is restricted to nearest- 
neighbour pairs. The site percolation corresponds to the 
q = 1 Potts model, and ferromagnetic Ising model to 
the q — 2 case. This model is one of the simplest non- 
trivial critical thermodynamic N-body system, and many 
exact or accurate results are known for standard values 
of (d,q). In particular, there exists a value q c (d) (for ex- 
ample : <7 C (2) = 4), for which when q < q c (d), then such 
an interacting system experiences a second-order phase 
transition at a finite critical temperature (for example : 
(3 C J = ln(l + y/q) at d = 2), while for q > q c {d) the 
transition is a first-order one. 

Let us consider here the case of the second order phase 
transition. All the scalings described above should hold. 



9 



We have first to define the order parameter for the sys- 
tem. If for a given configuration of the system, we call 
Nk the number of sites in the state k, where k varies from 
to q — 1, then the order parameter m is given by : 



q(N max /N) - 1 



1 



(54) 



Nmax in (|5j) is defined as the maximum of all Nk's. Fig. 
1 shows the m-distribution at the critical temperature 
in the (d, q) = (2,3) case, in the first-scaling form. The 
scaling is recovered very precisely, even for such small 
system sizes as 64 x 64. Note also the complicated shape 
of the scaling curve. 




(1) 

FIG. 1. Scaled m-distribution for the 3-states Potts model 
(q — 3) on 2-dimensional L x L square lattice, at the crit- 
ical temperature. Three sizes are shown : L = 32 (stars), 
L = 48 (circles) and L — 64 (squares). The thermalization is 
achieved after 4 x 10 6 Monte-Carlo steps, starting from the 
initial disordered system. 

We can discuss this scaling here in the slightly dif- 
ferent context of correlated variables jl7|]. Consider for 
simplicity, the Ising model (i.e. (d, q) = (d, 2) case). The 
extensive order parameter is just the sum of N correlated 
variables : M = s « ■ When the system is disordered, 
the spins are correlated on the short distance £ {£/N — > 
at the thermodynamic limit) and their mean value is zero. 
The central limit theorem tells then that the distribution 
of the random variable M/y/N is Gaussian when N be- 
comes large, with zero mean and finite variance. This 
can also be expressed by the asymptotic law : 



M 2 " 
'2 < M 2 >, 



(55) 



which is correctly in the first-scaling form, with a Gaus- 
sian shape and the trivial anomalous exponent g = 1/2 
(see (p|)). 



On the other hand, if the system is in the ordered 
phase, the average value of the individual spins is finite, 
say : < s, >= m, and the same reasoning can be used 
for the variable : (M — Nm)/yfN. This variable is of 
zero mean , finite variance and short-ranged correlated. 
So, its fluctuations are Gaussian and can be put in the 
second-scaling form : 



< M > 1/2 P[M] 



1 



exp 



(M— < M >f 
2 < M > 



(56) 



Of course, the most interesting case corresponds to the 
critical temperature. At this point, the spins are corre- 
lated throughout the system, and the magnetization can- 
not be evaluated by the central limit theorem. Instead, 
we can remark that the spin-correlations are power-law : 



1 



with r\ a critical exponent whose value should be between 
2 — d and 2. Looking at the total magnetization as the 
sum of N correlated variables, one gets 



< M 2 > = J2 < SiSj >= 



< s 2 > 



ho 



N J < s 5 Sf > r d ^dr 



(57) 



with L ~ N 1 /* 1 the typical macroscopic length of the 
system. This means that : 



< M 2 >~ 7V 1+ hr 



as the leading behaviour. This non-trivial anomalous ex- 
ponent 



n 



2d 



between 1/2 and 1, is here the sign of the criticality. 
The first-scaling law should hold in this case, as for the 
(d, q) = (2,3) Potts model discussed above, but the scal- 
ing function should be different since it depends on the 
precise form of the interactions. Only the tail can be 
linked to another critical exponent, as it has been writ- 
ten in Sect. II. C. 



V. REVERSIBLE AGGREGATION 
PERCOLATION MODEL 



THE 



Percolation model can be defined as follows. In a box 
(a part of the regular lattice), each site corresponds to a 
monomer and a proportion p of active bonds is set ran- 
domly between sites (the bond percolation model) . Such 
a network results in a distribution of clusters defined as 
ensemble of occupied sites connected by active bonds. 
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For a definite value of p, say p cr , a giant cluster almost 
surely spans the whole box. The sol-gel transition cor- 
responds to the appearance of 'infinite' cluster (gel) at 
a finite time. 'Infinite' in this context means that gel 
contains the finite fraction of total mass of the system. 
The sol - gel transition in finite systems can be suitably 
studied using moments of the number-size-distribution 
n s , i.e. the number of finite clusters of size s : 



M' k = 



S<Sma 



s k n s 



(58) 



where the summation is performed over all clusters with 
the exception of the largest cluster s = s max . The su- 
perscript ' recalls this constraint on summation in Eq. 
( p8| ) . The mass of the largest cluster is then : N — M[ , 
with : 



N 



E 



In infinite systems, one works with the normalized mo- 
ments of the concentration-size-distribution c s , i.e. the 
concentration of clusters of size s : 



E 



s k c. 



(59) 



where the summation in (|59|) runs over all finite clusters. 
Generally, concentrations are normalized such that : 

,. n s 
c s = lim — . 

N-*oo N 

The probability that a monomer belongs to the infinite 



cluster (gel) is equal to 1 



777 



1 ! 



with 



,. M 'k 
m k = lim — ^ . 

For example, in the thermodynamic limit when the size 
of the box becomes infinite , a finite fraction of the total 
number of vertices belongs to this cluster. Therefore, we 



= 1 for p < p cr and m'i < 1 for 
m! x is a decreasing function of the 



get the results : ml 
p > p cr . Moreover, 
occupation probability. This typical behaviour is com- 
monly (and incorrectly) called : 'the failure of mass con- 
servation 1 , but, as stated before, m\ is more simply the 
probability for a vertex to belong to some finite cluster. 



A. Percolation on the Bethe lattice 

The bond-percolation on the Bethe lattice with coor- 
dination number z, has been solved by Fisher and Essam 
|j39f . Here the main result we are interested in, is the 
concentration-sizc-distribution |39(1 : 



((z-l)s)\ 



2)!s! J 



_1 (1 - p)( 5 ~ 2 ) s+i 



and the first normalized moment 



l-p 



I-P* 

with p* being the smallest solution of equation : 

p*(l-p*f- 2 =p(l-pf- 2 . 
Let us define : 

1 



1 



For p < p cr , the only solution of the above equation is 
: p* = p, but when p is larger than p cri then there is a 
smaller non-trivial solution which behaves as p cr — \p— p cr \ 
nearp cr . Above this threshold, the moment m\ is smaller 
than 1 and behaves approximately as : 



1 



2(p-p cr ) 



1 



Pc 



The marginal case z — 2 corresponds to the linear-chain 
case. 

Coming back to the concentrations, we can see that for 
large values of the size s, the following Stirling approxi- 
mation holds : 



-5/2 



exp(— as) 



with a given by : 

a = In 



Pc 



1-Pc 



z-2 



For this model, a power-law behaviour of the concentra- 
tions is seen at the threshold p cr , namely : c s ~ s~ T with 
t = 5/2. Outside this threshold, an exponential cut-off 
is always present p(i| . This sort of critical behaviour 
at equilibrium is analogous to the thermal critical phe- 
nomena, and in particular, there exist two independent 
critical exponents, for example r and a. The latter one 
is the exponent of the mean cluster-size divergence. To- 
gether, those two critical exponents : r = 5/2 and er = 1, 
describe completely the critical features. 

This singular behaviour is due to the appearance of 
a giant cluster, the so-called percolation cluster, at the 
transition point. More precisely, in the infinite system 
the probability for a given site to belong to this infinite 
cluster is zero below the critical threshold p cr and positive 
above it. This probability is non- analytical at the critical 
point. Because of this behaviour, the extensive order 
parameter defined for finite systems is just the size of the 
largest cluster s max . 

As discussed in Sect. II. C, the corresponding finite-size 
order parameter scales as 



N 



2/3 
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Even though the system experiences the second-order 
critical phenomenon, fluctuations of the multiplicity dis- 
tribution remain small and the KNO-scaling does not 
hold. Of course, m' is not in this case an order parame- 
ter since r > 2 even though there is some irregularity in 
its behaviour passing the threshold. This non-analyticity 
can be illustrated by the exact result for bond-percolation 
on the Bethe lattice. In this mean-field case, the normal- 
ized O t/l -moment is : 



< = (1 - -p* 



2(z-l) 



l-p* 
(z-l)e 



2z-2 



with : e = p — p cri and e< 1. There is a jump of the first 
derivative of m' with respect to p : — z/2 for p — > p~ r , 
and (4 — 3z)/2 for p — > p+ . 



B. 3-dimensional percolation 

As shown by Botet et al [|3| , the multiplicity distribu- 
tion for the 3-dimensional bond percolation model on the 
cubic lattice at the infinite-network percolation threshold 
exhibits the A-scaling with A = 1/2 , and hence the frag- 
ment multiplicity is not related to the order parameter 
in this process. This is shown in Fig. 2c as a typical 
example of non-critical parameter scaling. Note that the 
multiplicity distributions in Fig. 2c are plotted in semi- 
logarithmic form to show clearly the Gaussian behaviour 
(parabolic shape on the figure). The proper order pa- 
rameter for this model is the normalized mass of the gel- 
phase, i. e. the mass of the largest cluster divided by the 
total mass of the system s max /N . Different probabil- 
ity distributions PN[s m ax/N] for different system sizes 
TV can be all compressed into a unique characteristic 
function (see Fig. 2. a) : 



< 



TV 



A' 



TV 



> 



< S r , 



> 



which is analogous to the KNO-scaling function. As an 
application of results developed in Sect. II. E, Fig. 2b 
shows the A-scaling for the shifted order parameter : 
M[ = TV — s max . The value of A (= 0.8), is consis- 
tent with the value of the anomalous dimension (25) : 
g = 0.8435, for the accepted values of the critical ex- 
ponents /3,7 in the 3-dimensional percolation p0[ . One 
should also remember, that A has been extracted from 
the small size (TV = 14 3 ,20 3 ,32 3 ) percolation network 
calculations at the infinite-network percolation thresh- 
old. This explains a small difference between the value 
for A from the scaling analysis and the expected value : 
A = g, in the infinite network. 

According to the results derived above for the second 
order phase transition, the second-scaling should hold 



outside of the critical point. This is correctly realized 
with the three variables s max , M[ and Mq for large or 
small values of the probability p. Fig. 3 shows such re- 
sults for the value p = 0.35. 
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FIG. 2. (a) The first-scaling of s m(la: -distributions at the 
percolation threshold (p = p cr = 0.2488) of the 3-dimensional 
bond percolation for lattices of different sizes : N = 14 3 (di- 
amonds), N = 20 3 (squares), TV = 32 3 (circles). The data 
corresponds to 10 5 events. 

(b) The A-scaling of the distributions of M[ = TV — s max for 
the same conditions as in (a). 

(c) The second-scaling of the multiplicity distributions plot- 
ted in log-linear scale (i.e. log(< Mq > 1,/2 P[Mo]) vs Z(i/2)) 
for the same conditions as in (a). 

Finally, it is instructive to see how the first-scaling is 
disappearing when the value of p is slightly shifted away 
from its critical value. Fig. 4 illustrates the deviations 
from the first-scaling for the values of parameter p close 
to p cr , on both sides of p cr . Even very close to the critical 
point, these deviations are quite significant and can be 
easily seen in this representation. 
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VI. IRREVERSIBLE AGGREGATION - 
EXAMPLE OF SMOLUCHOWSKI KINETIC 
MODEL 

The irreversible sol-gel transition can be modelled us- 
ing the coupled non-linear differential equations in dis- 
tributions c s of clusters of mass s per unit volume (the 
Smoluchowski equations |27 



dc s 
dt 



+j=s 



(60) 
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FIG. 3. (a) The second-scaling of s mal -distributions above 
the percolation threshold (p = 0.35) of the 3-dimensional 
bond percolation for lattices of different sizes : N = 14 3 (dia- 
monds), N = 20 3 (squares), N = 32 3 (circles). The calculated 
data corresponds to 10 events. 

(b) The second-scaling of the distributions of M{ = N — s max 
for the same conditions as in (a). 

(c) The second-scaling of the multiplicity distributions plot- 
ted in log-linear scale {i.e. log(< Mo > 1//2 P[Mo]) vs Z(i/2)) 
for the same conditions as in (a). 
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FIG. 4. (a) The s mal -distributions are plotted in the 
first-scaling form for parameters p close to the percola- 
tion threshold of the 3-dimensional bond percolation : (a) 
p = 0.252 ; (b) p = p cr = 0.2488 and (c) p = 0.245. The 
calculations are done for lattices of different sizes : TV = 14 3 



(diamonds), N = 32 3 
sponds to 10 J events. 



(circles). The calculated data corre- 



Coefficients K%j represent the probability of aggrega- 
tion per unit of time between two clusters of mass i and 
j . The Smoluchowski equations are derived from the 
master equation in the mean-field approximation |j3f : 

< c k ci >=< c fe >< q > . 

The time t includes both diffusion and reaction times. 
Eqs. ( J60| ) suppose irreversibility of the aggregation, i.e. 
the cluster fragmentation is excluded. One should notice 
however that the sum over j in (pOj) does not include the 
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infinite cluster (gel) because : 

c j=00 = l/oo 







Experimentally known aggregation kernels K^ are ho- 
mogeneous functions |28| : 

Kai,aj — ^i}3 

with A being the homogeneity index. Perhaps the sim- 
plest physically relevant example of a homogeneous ker- 
nel is : Kij — (ij)* 1 - It has been shown in this case 
that if fj, is larger than 1/2, then there exists a time t cr 
(t cr < oo) such that m! 1 = 1 for t < t cr but m\ < 1 for 

t --^ > t ct [P^IJ^m] ■ 

Let us consider now the case : Kij = [ij) 1 * with \l = 1 
in more details. It was shown in this case [^9| that the 
critical gelation time is : t cr = 1, and the solution for 
size-distribution of the Smoluchowski equations with the 
monodisperse initial condition is |32] | : 

s-2 

c, = -^e" 8 *^- 1 for t < 1 



s s 2 exp(— s) 



t 



for t > 1 



The asymptotic solution for large s is : 

c, l — S - b/2 exp [-s(t - 1 + In t)] for t<l 

~ 5/2 for £ > 1 (61) 



W2ir 
1 



Note that the power-law behaviour (r = 5/2) is present 
for t > 1 and not only at the threshold. The whole distri- 
bution of finite-size clusters evolves self-similarly and the 
appearance of a power-law behaviour is not here a sign 
of the critical behaviour but a specific characteristics of 
the gelation phase. 

The solution for the first normalized moment is : 



1 



t 



for t < 1 



for t>l 



With those asymptotic forms of c s one can calculate the 
gel fraction in the infinite system before and after the 
critical point : 



m G — 



m G = 1 



t 



for t < 1 



for t > 1 



It has been shown that gelation is analogous to a dy- 
namical critical phenomenon with [B7fl : 



lim — < S m ax > 



as the order parameter. For one realization, s m ax corre- 
sponds to the mass of the gel above t c = 1. 

For finite sizes, one makes the usual assumption that 
there exists a characteristic size diverging at the transi- 
tion, say : 



Nr. 



U-l 



-l/crjv 



(62) 



such that for the mass gel in the finite system one has 

< s max > ~ (t - 1)/ for t > 1 

In particular, at the gelation time one has : 



N g 



(63) 



Using the formula (29), which is valid both for the equi- 



librium and non-equilibrium systems, one can calculate 
the anomalous dimension. Given the value of r (see (|Sl[)), 
one finds g — 2/3. Hence : <jn — 1/3 can be deduced 
from (p3|). The average value of the order parameter 
< Smax > increases logarithmically for t < 1, and is a 
finite proportion of the system size when t > 1. 
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FIG. 5. The first-scaling of the s maa; variable in the Smolu- 
chowski kinetic model with the kernel Kij = ij at the critical 
time t = t cr — 1. The calculations are performed for two sys- 
tem sizes : N = 2 10 (diamonds) and N = 2 14 (circles). Each 
data set corresponds to 10 5 independent events. 

The illustration of the above discussion is shown in 
Figs. 5 and 6. Fig. 5 shows the distribution of s max in 
the first-scaling variables for systems of different sizes. 
The results have been obtained in the Smoluchowski 
model with the kernel : Kij = ij, at the critical time 
: t = t cr = 1- Fluctuation properties of s max outside 



of the critical time : t 



2t r 



are shown in Fig. 6. 
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The remaining parameters of the Smoluchowski calcula- 
tions are the same as used in the calculations shown in 
Fig. 5. In this case, the data for different system sizes 
collapse into the universal curve in the scaling variables 
with A = 1/2. One should keep in mind that in both 
cases, the fragment-size distribution is a power-law with 
t = 5/2 (sec @). 
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FIG. 6. The second-scaling of the s max variable in the 
Smoluchowski kinetic model with the kernel Kij = ij above 
the critical time t = 2t cr = 2. The calculations are performed 
for two system sizes : N = 2 10 (diamonds) and N = 2 14 (cir- 
cles). Each data set corresponds to 10 s independent events. 

logtfcCzj)] 








FIG. 7. The plot of the large - zm tail of the decimal loga- 
rithm of the scaling function : log^(z^), against zf^ for the 
system size N = 2 12 . The calculated data correspond to 10 5 
independent events. The solid line shows the dependence : 
3?(z(i) ~ exp(— Zni), which is expected from the value of the 
anomalous dimension g = 2/3. 

The relation between the form of tail of the scaling 



function and the anomalous dimension ( pq ) was derived 
analytically in Sect. II. C for the equilibrium systems at 
the second-order phase transition. For non-equilibrium 
systems, we do not know equally rigorous derivation (see 
also Sect. VI. A). On the other hand, one may expect that 
the relation between the TV-dependence of the average 
value of the order parameter and the asymptotic form of 
the scaling function in the limit N — > oo, i.e. between v 
and g, is connected to the asymptotic stability of the limit 
distributions. Actually, there is a very close connection of 
the renormalization group ideas and the limit theorems 
in the probability theory j34j|. If true, then the relation 
( pfj| ) could be valid in a more general framework than the 
one provided by the equilibrium statistical mechanics. 
To check this assertion, in Fig. 7 we show the plot of 
the logarithm of the scaling function (see Fig. 5) 

versus for large values of Z(iy If the relation ( |26| ) is 
valid also for the non-equilibrium sol - gel second-order 
phase transition, then : ~ exp(— z 3 ), and the tail 

of the scaling function should be a straight line in Fig. 
7. That is indeed the case. 

Figs. 8 and 9, show the A-scaling for the shifted order 
parameter variable : Mj = N — s max . Results of the 
Smoluchowski calculations with the kernel : Kij = ij, 
are shown at t = t cr (see Fig. 8), and at t = 2t cr (see 
Fig. 9). One sees that the M x distribution exhibits a 
qualitative change while going from the critical time t = 
t cr , where A = 0.67, to t = 2t cr for which A = 1/2. At 
t = t cr , the value of A obtained by superposing different 
M 1 distributions in the scaling plot ( pf^ ) , agrees perfectly 
with the value of the anomalous dimension g (=2/3). 



<£(Z(0.67)) 




FIG. 8. The A-scaling of the distributions of M{ shifted 
order parameter in the Smoluchowski kinetic model with ker- 
nel Kij — ij at the critical time t = t cr — 1. Two system sizes 
are considered : iV = 2 10 (diamonds) and N = 2 14 (circles). 
The calculated data correspond to 10 5 independent events. 

By comparing Figs. 5, 6 with Figs. 8, 9, one may see 
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also that the effect of changing the variable : s max — > 
M l5 is seen only at the critical time (compare Figs. 5 
and 8) where (A = 1) — > (A = 0.67), and is absent 
above the critical time (compare Figs. 6 and 9) where A 
(=1/2) remains unchanged. 




^(1/2) 

FIG. 9. The second-scaling of the distributions of M[ 
shifted order parameter in the Smoluchowski kinetic model 



with kernel Ki, 



above the critical time, at t 

-,10 



Two system sizes are considered : N = 2 1U (diamonds) and 
TV = 2 14 (circles). The calculated data correspond to 10 5 
independent events. 
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FIG. 10. The distributions of M[ shifted order parameter 
in the Smoluchowski kinetic model with kernel Kij — ij at 
the critical time t = t cr = 1 , are plotted in the scaling 
variables of second-scaling. Two system sizes are considered 
: N = 2 10 (diamonds) and N = 2 14 (circles). The calculated 
data correspond to 10 s independent events. The failure of the 
second- scaling is clearly visible. 

Finally, in Fig. 10 we show the size-dependence of 



the Maldistributions at t = t cr , when the distributions 
are plotted in the 'wrong' variables of the second-scaling 
A = 1/2. The distributions for two system sizes are 
clearly displaced , showing the sensitivity of the scaling 
analysis and failure of the second-scaling. 



A. The origin of fluctuations and the argument of 
Van Kampen 

Q - expansion is a systematic expansion of the master 
equations in powers of 1/N ||. Lushnikov was the first 
to express the generating functions as the contour inte- 
grals for quantities like the moments M' k |35| . Then Van 
Dongen and Ernst j3£| used £1 - expansion to calculate 
explicitly these integrals for the moments M' k in some 
simple cases like : Kij = ij. For example, the result for 
M[ can be expressed in terms of the generating function 
for the s maar distribution PN[s ma x], as : 



5> 



N\S r , 



2iTr' 



x exp 



-Nu 



■ N 

E 



dz 

7n+i 



,s(N-s)/(2N) 



N 



s-l 



Using then the identity : 

gn 




Enl 
l! a iai!. ..n! a »a„ 



exp 



, where the sum runs over different sets {ai, ...a„} with 
the constraint : a± + ... + a n = n, and the particular 
result written down by Van Dongen and Ernst for the 



Ki 



zj-case 



exp 



" N 

E 

.8=1 

N 



N 



s-l 



exp[ — — - — \z 



2N 



EZ' 11 
-exp[--st- 

s— 1 

we can find : 

^ PN[s m ax] exp (S max u) 



X 



0{l/N) 



exp 



* E 



\k=l 



m' k (-u) h 



(64) 



Having the f2 - expansion of the generating function of 
the s maa; -distribution, we can conclude about the scal- 
ing at the gelation point. The moments m' k of the size- 
distribution for infinite systems are known to diverge near 
the gelation time p(| as : 
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m' k -\t-ir 2k . 
For finite systems, using : er/v = 1/3, one obtains : 

m' k ~\t-l\- 2k + s f k (N\t-l\ 3 )~N 2J ^ (65) 

at the gelation time. We have then found the asymptotic 
result : 



Nm' k ~ a k < s T 



where a k J s are some positive constants. Inserting ( p6[ ) in 
( |64| ) , one can show that the generating function for s max 
is the function of a single variable : < s max > u, what is 
a sufficient condition for the validity of the first-scaling 
law (g). 

We can have informations on similar scalings for var- 
ious moment-distributions. Q - expansion leads to the 
results : 



< My > 



< ML > 2 



Nm' 2k - (1 - t)m' fc 2 +1 



< M' k > = Nm' k 



(66) 



for the values of k when all the quantities are defined. At 
the transition (t = 1), the relation ( |65| ) allows to calculate 
the reduced moments m' k . The results can be written in 
the compact form : 



< M' k 2 > — < M' k > 2 
< M' k > 2A 



const 



with the following values of exponent A : 

A = 1/2 for k < 3/4 

A = 2fc/3 for 3/4 < k < 3/2 

A = l for 3/2 <k 

These are indications for A-scaling according to the re- 
mark ( f49| ) in Sect. II. H. More precisely : the moments of 
order k < 3/4 are not critical (the second-scaling law), 
the moments of order k between 3/4 and 3/2 exhibit the 
A-scaling with A = 2fc/3. In particular, for k = 1, one 
recovers the correct value : A = g = 2/3, corresponding 
to the general argument of the shifted order-parameter 
@ with a% = a% = 1. Finally, when the value of k is 
larger than 3/2, we obtain the first-scaling law for the 
distribution of moments M' k . This is also a consequence 
of the shifted order parameter argument, since in these 
cases : 

<M' k > ~ <s k rnax > . 

Far from the critical point, all the reduced moments m' k 
are independent of N , since the correlation size : 

1 



(see (pl|)) is finite. Then, for any value of k, the second- 
scaling law holds, as expected from the general theory. 

The above results about A-scaling for various moments 
of the size-distribution in the Smoluchowski model with 
kernel Kij = ij, are not complete since the arguments in- 
volve only the second cumulant moment K2- In principle, 
as shown in Sect. II. H, all cumulants should be inves- 
tigated. So, even though many exact results are known 
in this model, the complete analytical solution is not yet 
available. 

The same study as presented above for gelling systems, 
can be performed also for non-gelling systems. An ex- 



ample of this kind is obtained for : Ki 



j. In this 



case, the size-distribution is power-law with the exponent 
t = 3/2 |Q and, following the discussion in Sect. II. C, 
the cluster multiplicity can be the order parameter. One 
can derive analytically, that the multiplicity distribution 
is binomial : 



P N [M ,t 



(1 - exp( 



N — 1 
M Q -I 

x exp[-(M — l)Nt] 



-Nt)f- Mo x 



and can be approximated for N oo and for a finite 
value of < M > /N by : 



<M > 1/2 P N [M ,t] 



1 



2tt(1 - e- m ) 



x exp 



1 



(M— < M >Y 



2(1 - e-^*) < M > 



(67) 



t-1 + lnt 



what corresponds to the second-scaling. One may no- 
tice, that this binomial distribution is exactly equivalent 
to the bond percolation on a Bethe lattice with the oc- 
cupation probability : 

p = 1 - cxp(-Nt) . 

In spite of self-similar features in the fragment-size distri- 
bution at the infinite time, one does not see any critical 
behaviour in the cluster multiplicity distribution at any 
time in the non-gelling aggregation systems. This con- 
firms the observation made in Sect. III. A in the Mekjian 
equilibrium model, that the power-law size-distribution 
alone does not guarantee that the system exhibits the 
critical behaviour. 

The insight gained from the numerical simulations of 
Smoluchowski equations and the evidences from exact 
results for both gelling and non-gelling aggregation sys- 
tems, provide strong hints that the discussion of Sect. 
II. H is valid not only for the equilibrium systems but 
also for the non-equilibrium ones. We see the same sig- 
nificance of the A-scaling in non-equilibrium systems as 
found in thermodynamic systems, and not only at the 
critical point but also close to the critical point or even 
far from it. We believe that this universality, which is 
common to equilibrium and non-equilibrium systems, has 
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deeper foundation in the relation between renormaliza- 
tion group ideas for self-similar systems and the limit 
theorems of probability theory for the asymptotic scal- 
ing laws of order-parameter distributions. The concept 
of statistical equilibrium does not intervene at this level. 
One should also remember that the universality discussed 
in this work, is associated with only one critical exponent, 
and certainly does not exhausts all the singularity prop- 
erties of the thermodynamical potential in second-order 
thermal phase-transition. 

VII. CONCLUSIONS 

We have presented in this paper the theory of univer- 
sal scaling laws of the order parameter fluctuations in 
any system in which the second-order critical behaviour 
exist. These scaling laws, called A-scaling laws, are rig- 
orously derived for the equilibrium systems. Moreover, 
both analytical and numerical evidence is presented in 
favour of the general validity of the A-scaling laws also 
for the off-equilibrium processes which exhibit the critical 
phenomenon of the second-order. 

In this work we have discussed different aggregation 
models, both reversible and irreversible, finding the same 
connexion between scaling function properties and the 
anomalous dimension (critical exponents). These results 
can be important in the phenomenological analysis of 
'critical behaviour' in finite systems, where the critical 
exponent analysis is dubious and, moreover, the precise 
mechanism of the process may be unknown. In these 
cases, the A-scaling analysis allows to select both the 
relevant observable and the interesting initial conditions, 
which lead to the 'pseudocritical' behaviour in the stud- 
ied process. Another interesting aspect of the A-scaling 
analysis is the possibility of compressing data and, hence, 
the elimination of redundant dependences in the data on 
parameters like the system size, the total energy, or the 
total momentum etc.. This provides an obligatory in- 
termediate step in any phenomenological analysis before 
the laws governing complicated dynamics can be found. 
Finally, one should stress that the scaling laws discussed 
in this paper are independent of whether one deals with 
an equilibrium or an off-equilibrium process. This is a 
crucial advantage in the studies of short-lived systems. 
Examples of the multifragmentation processes in the col- 
lisions of atomic nuclei or atomic clusters well illustrate 
this problematic In the absence of thermal equilib- 
rium, which is a theoretical hypothesis difficult to verify 
in dynamically formed short-lived systems, we simply do 
not have at our disposal any other tool to address reli- 
ably the question of possible 'criticality' of the studied 
process. 

As said before, the A-scaling analysis developed in 
this work provides an alternative to the critical expo- 
nents analysis in the equilibrium systems and is the only 



tool for the analysis of the non-equilibrium systems. All 
essential information can be deduced from the scaling 
function, the value of A parameter, the form of the tail 
of the scaling function and the value of the anomalous 
exponent. With these informations it is possible to find 
out whether the studied process is at the critical point, 
in its neighbourhood or far away from it. Reference point 
in this analysis is the (approximate) self-similarity of the 
system. Generalization of the above scaling theory to 
discontinuous phase-order transitions for which the char- 
acteristic length can be defined, is in progress. 

We thank R. Paredes V. for providing us with the Potts 
model data which were used in preparing Fig. 1. 
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